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Abstract 



The Coulomb deexcitation of light exotic atoms in collisions with hydrogen atoms has 
been studied in the framework of the fully quantum-mechanical close-coupling method for 
I the first time. The calculations of the Z-averaged cross sections have been performed for 

{^p)n and {^id)n atoms in the states with the principal quantum number n = 3 8 and 
^SJ ■ relative energies region E = 0.01 100 eV. The obtained results reveal the new n and E 

^ ! dependences of the Coulomb deexcitation cross sections. The large fraction (up to ~ 36%) 

\ of the transition with An > 1 is also predicted. 

o ' 

^ 1 Introduction 
o 

In this paper we present ab initio quantum-mechanical treatment of non-reactive scattering pro- 
cesses of the excited exotic hydrogen-like atoms: 

^ \ - elastic scattering 

S^. ; {ax)n + {he')u {ax)n + {he')u] (1) 



■ " Stark mixing 

H : {a'X)ni + {he')v ^ {ax)nv + {he')u\ (2) 

- Coulomb deexcitation 

{ax)ni + {he')u {ax)n>u + {h'e)^. (3) 



Here (a, 6) = {p,d,t) are hydrogen isotopes and x = fi^ , , K~ , p; {n,l) are the principal and 
orbital quantum numbers of exotic atom and u are the hydrogen atom quantum numbers. The 
processes (1) - (2) decelerate while Coulomb deexcitation (3) accelerates the exotic atoms, influ- 
encing their quantum number and energy distributions during the cascade. The last process has 
attracted particular attention especially after the "hot" up atoms with the kinetic energy up to 
200 eV were found experimentally [1,2]. 

Starting from the classical paper by Leon and Bethe [3] , Stark transitions has been treated in 
the semiclassical straight-line-trajectory approximation (see [4] and references therein). The fully 
quantum-mechanical treatment of the processes (1) - (2) based on the adiabatic description was 
given in [5-8]. Recently [9,10], the elastic scattering and Stark transitions have also been studied 
in a close-coupling approach with the electron screening effect taken into account by the model. 
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Thus, at the present time we have more or less reahstic description of the processes (l)-(2) in a 
wide range of n and E. 

Concerning the acceleration process (3), the situation here is much less defined, especially 
for low n. The first work (to the best of our knowledge) on the CD process was performed by 
Bracci and Fiorentini [11] using the semiclassical approach with some additional approximations 
concerning the assumption of the large values of the quantum numbers under consideration (n ^ 
3) and some others. In the following numerous papers [12, 13] (and references therein) the CD 
process is considered within the asymptotic approaches using the adiabatic hidden crossing theory 
[14] for the pure Coulomb three-body systems or for the screened Coulomb centers chosen in the 
static form (the dipole screening). Recently, the calculations of the process (3) were also performed 
in the classical-trajectory Monte-Carlo (CTMC) [15] approach. While the Coulomb deexcitation 
cross sections obtained in CTMC approach are in fair agreement with the semiclassical ones of 
Bracci and Fiorentini [11], the much more elaborated advanced adiabatic approach (AAA) [12,13] 
gives much smaller Coulomb deexcitation cross sections which can not explain the experimental 
data. The reasons of such a strong discrepancy are not clear. One can only assume that the 
semiclassical model [11] as well as the CTMC approach are not valid for low- lying states. Thus, 
until now there is no satisfactory description of this process in the most interesting region (n = 
3^7). 

The main aim of this paper is to obtain the cross sections of Coulomb deexcitation for 3 < n < 8 
in realistic quantum-mechanical approach which is free from the additional approximations used in 
previous studies, in particular, the two-state approximation. For this reason we use a unified treat- 
ment of elastic scattering. Stark transitions and Coulomb deexcitation within the close-coupling 
method. This approach was used previously for the study of the elastic scattering and Stark tran- 
sitions of the excited muonic hydrogen in the collisions with the hydrogen ones [16] and a good 
agreement as compared with the results of the quantum-mechanical adiabatic description [5-8] was 
obtained. In the following section wc briefly describe the close-coupling formalism and the main 
expressions. The results of the closc-couphng calculations concerning the total cross sections of 
the process (3) are presented and discussed in Sec. III. Finally, summary and concluding remarks 
are given in Sec. IV. 



2 Theory 

The non-relativistic Hamiltonian for the four-body system {afi^ + be~), after separating the center 
of mass motion, can be written in Jacobi coordinates (R, p, r) as 

H = -^Ar + hM + he{r) + V{r, p, R), (4) 

where m is the reduced mass of the system, V is the electrostatic interaction between the subsys- 
tems, hf^ and hg are the Hamiltonian of the free exotic and hydrogen atoms 

K^nlmip) = £n^nlm{p): (5) 

he<Pis{r) = eis(pis{r), (6) 

where ^nim{p) ^-nd (pis{r) are the wave functions of the exotic atom and hydrogen atom bound 
states. En and ei^ are the corresponding eigenvalues. The interaction potential of the subsystems 
is defined by 

V{r, p, R) = Vab + V^, + Vae + V^e (7) 
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with the two-body Coulomb interactions 



Vab^ — ^\R + up- Uer\ \ V/^b = — — = -|R - - J^er\ \ (8) 

Tab 

V^, = — = \R-^p + ^,v\-\ Ke = -— = -|R + i/p + ^,v\-\ (9) 

Here we use the set of Jacobi coordinates (R, p, r): 

R^R^-R^a, p = r^-r„, r = re-rft, 

where r^, rf,, r^, rg are the radius-vectors of the nuclei, muon and electron in the lab-system , 
R//, R^a are the center of mass radius- vectors of the hydrogen and exotic atoms, respectively. 
The coefficients i/, ^, and in the two-body interactions (8)- (9) depend on the masses of the 
particles, 

v = m^l {nif, + m„), i = nia/ {rrif, + nia), (10) 
Ue = me/{me + fUb), Ce = rUb/ {m^ + rUb), (11) 

and satisfy 

i^ + e = ^e + ee = l- (12) 

{■ma,mb,mn and me are the masses of the hydrogen isotopes, muon and electron, respectively). 
Atomic units (a.u) h — e — meinb/ime + mb) — 1 will be used throughout the paper unless other- 
wise stated. 

The total wave function ^(p, r,R) of the system satisfies the time independent Schrodinger 
equation with the Hamiltonian (4): 

(-^Ar + h^ + he + V)-^{p, r, R) = E-^{p, r, R), (13) 

where E is the total energy of the system. 

In this paper, as well as in the previous studies [11, 12, 15], we assume that the state of the 
target electron is fixed during the coUision. The electron excitations can be taken into account 
in a straightforward manner. Owing to the rotation and inversion symmetries of the total sys- 
tem, the vector solutions of (13) is introduced the total angular momentum representation. In 
a space-fixed coordinate frame we built the basis states from the eigenvectors of the operators 
he, hfj,, 1^, L^, J^, Jz and the total parity tt with eigenvalues Su, En, l{l + 1), L{L + 1), J( J -|- 1), M 
and (—1)'+^, respectively: 

|r)^0i.(r)|7) (14) 

where 

|7) = \nl, L:JM)= ^(ZmLA| JM)$„^^(p)yi,A(R), (15) 

\T) = -^R,,{r)Rni(p)y^''{p,R), (16) 
v47r 



yi^'ip, R) = J2(^mLX\JM)YUp)YL^{R) (17) 



mX 
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Here the orbital angular momentum 1 of {aiJ,)ni is coupled with the orbital momentum L of the 
relative motion to give the total angular momentum, J = 1 + L. Then, for the fixed values of 
J, M, TT = (— 1)'+-^ the exact solution of the Schrodinger equation for the coUiding system, 

(E-i7)^^^-(r,p,R)=0, (18) 

is expanded as follows 

vl/r^(r,p,R) = ^u{r)^J2'^nfL{R)\nl,L: JM), (19) 

nlL 

where the G:ljj^{R) are the radial channel functions and the sum is restricted to (/, L) values 
to satisfy the total parity conservation. This expansion leads to the coupled radial scattering 
equations 

+ " ^^^^W^) ^^'^^^^ = 2^ E ^nlL,n'l'L'iR) G'J.M, (20) 
^ ^ n'l'L' 



where k"^ = 2m{E — — eis) = 2m{Ecm + ^£nn') specify the channel wave numbers; Ecm is the 
relative motion energy in the entrance channel, A£„„/ = 0.5/x(n^ — {n'Y)/{nn'Y is the difference 
of the exotic atom bound energies in the initial and the final states, // is the reduced mass of the 
exotic atom. Finally, W^u^ ^ntj^, are the matrix elements that couple asymptotic channels {nlL] J) 
and {n'l'L'; jy. 

"^Ll^u'W = (ls,n/,L : J\V\\s,v!l' ,U : J). (21) 
The radial functions G'^E,n'VL'iR^ must be regular everywhere, and, at — > 

Gtn'VL'{^) = 0(~ it:^+^) (22) 
and at asymptotic distances {R — > oo) satisfy the usual boundary conditions 

Gi^n'L'iR) =^ -^{SifSnn'Su'SLue-'^'-''-''''/'^ - S'inl, L ^ n'l', L')e^(>^fR-L'-/^)}^ (23) 



where ki, kf are the wave numbers of initial and final channels and S'^{nl,L — > n'l'.L') is the 
scattering matrix in the total angular momentum representation. The indexes of the entrance 
channel and target electron state are omitted for brevity. The scattering amplitude for nlm — > 
n'l'm' is defined by 

f{nlmL n'l'm' L'\ki,kf,R) = V i^'-^{lmLX\JM){l'm'L'X'\JM)x 



ki kf 



JMLL'XX' 

J I 



X Yl^{]<i^YL^y{K)T\nlL ^ n'l'L'), (24) 
where the transition matrix T"^ used here is given by 

T'{nl, L ^ n'l', L') = 5nn'h>hL>5mm'hx> - S'{nl, L ^ n'l', L'). (25) 
The matrix elements (21) of the interaction potential (7-9) are given by 

WnlL,n'l'L'{R) = ^J dr dp dRRUr)Rnl{p)Rn' I' (p) 

X (3^r )*(P, R) V{r,p, R) yi^^ip, R), (26) 



4 



where the radial hydrogen-hke wave functions are given exphcitly by 

I n-l-l 



Rnlip) = N^i exp(-p/na) ^ S,{n,l) (J^J (27) 



q=0 

(a is the Bohr' radius of the exotic atom in a.u.) with 



\na J 



{n + iy.{n - I - ly. 

2n 



1/2 



(28) 



^>'^)^H- ^!(n-/-l-'.)!(2/ + l + ,)r ^2^) 
Averaging V{r, p, R) in (26) over Is-state of hydrogen atom we obtain 

1 /"^ 

V{R,p) = -j^ drRUr)V{r,p,R)^ 

= hu.^^^{R,p) - [/-s,«.(R,p)} - -^{t/.,..(R,p) - ^-c,«.(R,p)}, (30) 



where 



2x|R+ap| 



3 2|R+ap| / 1 9 \ e " 

Ua0('R,p) = (1 + r)e"^^ = limll---^] pj -. (31) 



Using the additional theorem for the spherical Bessel functions 

„-A|Ri+ri| A 



Ri + ri 



X 



AVi/2(Ai?i) /t+i/2(Ari) \r,<R, + It+i/2{XRi) i^i+i/2(Ari)|^^^^J (32) 



{Ip{x) and K.p{x) are the modified spherical Bessel functions of the first and third kind) and 
substituting eqs. (30-32) into (26), we integrate over the angular variables. Furthermore, using the 
angular momentum algebra and integrating over p, we obtain: 



= (-l)-^^'+''/-^^'-'-^v1l>IZ^£(/0/'0|tO)(LOL'0|iO) ^ 



I I' t 



t=0 



\ L' L J 



X \ i [{-lfWt{R,v,i,-M,n'l')-Wt{R,i,ie-M,n'l')] 



- [{-rfWt{R, V, ve\ nl, n'l') - Wt{R, i^e] nl, n'l')] \ (33) 



{tm is the maximum value of the allowed multipoles). Here the next definitions are used: 

Wt{R, a, (3; nl, n'l') = Ki,n'i' V S^, (n, V S^, {n', I') x 

X {Ht{z)Ji''{z, A(n, n', a, /?)) - ht{z)J2''{z, A(n, n', a, /5)) + 
+ Ft{z)f/{z, A(n, n', a, /?)) +ft{z)f^\z, A(n, n', a, /?))} , (34) 
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where z = 2R//3, s = I + I' + rui + 1712, L = 2L + 1; 

'+1 / o„ \ ''+1 



Ki,n'i'^^^(^^] [^^] ^{n + l)\{n-l-l)\{n' + l')\{n'-l'-l); (35) 
n + n'\n + n'J \n + n'J 



. , . 2nn' aa 
X{n,n,a,p)^ — ■ — - — ; (36) 
n + n p 



and 



Ht{x) = (1 - 2t)ht{x) + xht+iix); (37) 

Ft{x) = (1 - 2t)Mx) - xft+,{x). (38) 
The functions ht{y) and ft{y) are given by 

ht{y) = \f^Kt+i/2{y) (39) 
V Try 

ft{y) ^ ^lYy^^^'/'^y)- (40) 

The radial integrals J*''*(a;, A) are defined as follows: 

Jl'^(x, A) = y'^^e-yft{\y) d|/, (41) 
^0 

J*'^(x,A) = AJl+^'^+^(x,A), (42) 

J*'^(a;,A)= / y'^\-yht{\y)dy, (43) 

J*'^(a;,A) = AJ*+^'^+\a;,A) (44) 

and calculated analytically using the power series for the modified Bessel functions. 

Finally, we give the exphcit expressions for the cross sections to be discussed in this paper. 
The partial cross sections of the processes (1) - (3) for the transitions nl — > n'V , averaging over an 
initial distribution of the degenerated substates and summed over the degenerate final substates, 
are given by 

<j\nl n'l'; E) = Yl l^'H^ ^ ^'^'^Ol'- (45) 

The /-averaged partial cross sections for the transitions n ^ n' are then computed by summing 
over I and I' with the statistic factor (2/ + l)/n^: 

'^Ln'iE) = ^^-^ E \T'inlL - n'l'ir. (46) 

* 1,1' LL' 

The total cross sections for the transition n — > n' is obtained by summing the corresponding 
partial cross section over the total angular momentum J: 

ann'(E) = Y,^L'(E). (47) 
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For the discussion of the obtained results and the comparison of them with the other approaches 
we will also need the total cross section of Coulomb deexcitation including all transitions n — > n' 
with n' < n: 

<J^''{E)^Y.^nn'{E). (48) 
n'<n 

3 Results 

The close-coupling method described in the previous Section has been used to obtain the cross 
sections for the collisions of the and ji'd atoms in excited states with hydrogen atoms. The 
present paper had at least two goals: first, to apply the fully quantum-mechanical approach for 
the study of the processes (1) - (3) and, second, to clear the validity of the energy and principle 
quantum number dependence of the Coulomb deexcitation cross sections used in literature and 
based on the semiclassical model [11]. So, we present here only a small part of our results. In 
particular, to illustrate some main features of the calculated cross sections, we discuss only the 
/-averaged CD cross sections. The detailed results of the calculations will be published elsewhere. 

The coupled differential equations (20) are solved numerically using the Numerov method and 
with the real boundary conditions involving the X-matrix instead of jS-matrix. The corresponding 
T-matrix can be calculated from equation 

T = 2K{I - iK)-^ = 2X^(1 + K'^)-^ - 2iK{l + X^)"\ 

where / is the unit matrix. In the calculations we have taken into account all possible multipoles 
in the interaction potentials and all open channels with n' < n. 

The close-coupling calculations have been carried out for the relative collision energies 
from 0.01 up to 100 eV and for the excited states with n = 3 -r 8. At all energies we take into 
account so many values of the total angular momentum J and the relative angular momentum L 
that all the inelastic cross sections were calculated with the accuracy not less than 0.1%. 

Table 1: Dependence of the elastic and deexcitation cross sections cr^^/ for (/i(i)„ + H coUisions 
at energy Ecm = 0.1 eV on the number tm of terms included in the multipole expansion for the 
matrix elements (33). 





<^66 


0"65 


0"64 


0-63 


0-62 


(761 


1 


54.70 


2.36 


0.20 


0.012 


0.0003 


0.0000002 


2 


50.24 


5.55 


0.37 


0.055 


0.0052 


0.000008 


3 


51.58 


4.82 


0.38 


0.038 


0.0176 


0.000038 


4 


52.65 


2.71 


0.38 


0.031 


0.0177 


0.000060 


5 


52.31 


2.72 


0.39 


0.026 


0.0277 


0.000079 


6 


54.19 


1.42 


0.36 


0.031 


0.0295 


0.000093 


7 


54.77 


1.04 


0.32 


0.032 


0.0279 


0.000100 


8 


54.46 


1.18 


0.30 


0.030 


0.0270 


0.000102 


9 


54.36 


1.17 


0.29 


0.028 


0.0269 


0.000103 


10 


54.31 


1.20 


0.30 


0.029 


0.0271 


0.000103 



The results of the calculations arc presented in Tables 1-3 and Figs. 1-6. The effects of the 
higher multipoles in the expansion of the interaction potentials are large (see Table 1). In contrast 
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to the elastic scattering and Stark transitions in which the "dipole" approximation (tm = 1) gives 
practically the exact results, there is no case for the Coulomb deexcitation. Ignoring the higher 
multipoles would lead to the significant distortions of the cross sections. As it is seen from 
Table 1 the results are sharply varied in the top part of the table (t^n = 1^7) and are practically 
unchanged in the bottom one (tm = 8^10). Thus, to provide the proper treatment of the Coulomb 
deexcitation all the multipoles in the expansion of the interaction potential must be taken into 
account. 

In Fig. 1 we present the J dependence of the partial- wave cross sections a^^, ior n — 7 at three 
fixed energies 0.1, 2, and 50 eV. It is seen that a substantial part of the Coulomb cross sections 
comes from the partial waves with rather a low J (in comparison with the elastic and Stark mixing 
processes). It is worthwhile to note that for n = 3 the range of J values is much less, e.g. at 
relative energies 0.01, 2 and 100 eV we obtain = 7, 11, 15 respectively, whereas for the elastic 
scattering we obtain J^ax = 7, 60, 120 for the same energies. 




J 

Figure 1: The partial cross sections of the Coulomb deexcitation (T^ g and the elastic scattering o"y 7 
in a.u. for (/xp)„=7 + H collisions as function of the total angular momentum J at three energies: 
0.1 eV (dashed), 2 eV (dotted) and 50 eV (solid). 

The convergence of the results upon the number of open channels (the transitions with An > 1 
are included) is illustrated in Table 2 for {fj,p)n=7 + D collisions at the relative energy Ecm = 1 eV. 
It is seen that only the nearest channels are most important. While n increasing the energy 
gap between the exotic atom states decreases. As a result the infiuence of the nearest channels 
becomes more significant and leads to the new n and E dependence of the Coulomb deexcitation 
cross sections as compared with the semiclassical ones. 

The energy dependence of the total Coulomb deexcitation cross sections a^^{E) for n = 3 -^ 8 
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Table 2: The dependence of elastic and Coulomb deexcitation cross sections cr7„/ for (/ip)„ + D 
collisions at energy Ecm = 1 eV on the number of open channels iVievei included in the calculation 
(all terms of the multipole expansion for the potential are taken into account). 



^Icvel 


(T77 




(^75 


0-74 


(^73 


(^72 


0-71 


2 


39.578016 


1.063412 












3 


40.181324 


0.679067 


0.096652 










4 


40.181170 


0.598721 


0.056935 


0.009409 








5 


40.177221 


0.598219 


0.055336 


0.009368 


0.001601 






6 


40.177061 


0.597370 


0.055110 


0.009343 


0.001635 


0.001241 




7 


40.177065 


0.597360 


0.055110 


0.009343 


0.001635 


0.001241 


0.000005 



is shown in Fig. 2 in comparison with the results of Bracci and Fiorentini [11] for n = 3, 5, 7 and 
the CTMC calculations of Jensen and Markushin [15] for n = 4 8 and Ecm = 0.5 and 5 eV. The 
present results are in poor agreement with the semiclassical results [11] and CTMC results [15]. 
There is only a qualitative agreement for n = 5 and energies above 2 eV. The energy dependence 
of the Coulomb deexcitation cross sections in our approach is approximately given by 1/ \fE^ at 
-E'en! 0.5 eV (see Fig. 3), in contrast to the IjE^m behavior found for low energies in [11] and also 
in the advanced adiabatic approach (AAA) [13]. 




0.01 



0.01 



0.01 



0.1 



10 



, eV 



Figure 2: The cross sections of Coulomb deexcitation cr^^ in a.u. for (fip)n+H collisions calculated 
in the quantum-mechanical close-coupling method (solid lines) in comparison with the results of 
Bracci and Fiorentini [11] (dashed) and classic Monte-Carlo results [15] (triangles). 



As for the n dependence of the total Coulomb deexcitation cross sections a^^{E) the situation 
is more complicated than it is followed from the semiclassical model [11] (see Fig. 2). There is no 



dependence 



CD 



n ' 



(7 ~ 2) as it is followed from [11]. Our prediction for the n dependence 



of the CD cross sections is quite different. In our opinion, this behaviour is explained by the 
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Table 3: The cross sections of Coulomb deexcitation for {fip)n + H collisions calculated in the 



quantum-mechanical close-coupling method. 



Ecm, eV 


0.01 


0.05 


0.1 


0.2 


0.5 


1.0 


2.0 


5.0 


7.0 


10.0 


15.0 


20.0 


(^3,2 


4.821 


1.505 


1.032 


0.803 


0.403 


0.209 


0.105 


0.043 


0.032 


0.023 


0.015 


0.012 


CD 


4.848 


1.511 


1.035 


0.805 


0.403 


0.210 


0.105 


0.043 


0.032 


0.023 


0.015 


0.012 


/a, % 


0.6 


0.4 


0.3 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


(^4,3 


4.126 


1.730 


0.936 


0.704 


0.515 


0.333 


0.189 


0.083 


0.060 


0.043 


0.030 


0.024 




5.526 


2.656 


1.397 


0.965 


0.637 


0.403 


0.226 


0.098 


0.071 


0.052 


0.036 


0.029 


/4, % 


25.3 


34.9 


33.0 


27.0 


19.1 


17.2 


16.4 


16.0 


16.2 


16.3 


16.5 


16.7 


0"5,4 


4.364 


1.349 


0.943 


0.574 


0.306 


0.211 


0.147 


0.070 


0.054 


0.041 


0.031 


0.026 




4.955 


1.705 


1.270 


0.819 


0.422 


0.277 


0.182 


0.085 


0.065 


0.049 


0.037 


0.031 


/5, % 


11.9 


20.9 


25.8 


30.0 


27.4 


23.9 


19.0 


17.2 


16.5 


16.1 


15.4 


14.9 




4.063 


1.804 


1.619 


0.856 


0.442 


0.250 


0.136 


0.058 


0.043 


0.032 


0.024 


0.021 


CD 


6.246 


2.598 


2.166 


1.189 


0.661 


0.387 


0.211 


0.090 


0.067 


0.050 


0.038 


0.033 


/6, % 


35.0 


30.6 


25.3 


28.0 


33.1 


35.4 


35.5 


36.0 


36.5 


36.8 


36.8 


36.1 




5.954 


2.888 


2.005 


1.410 


0.902 


0.488 


0.275 


0.119 


0.088 


0.065 


0.048 


0.040 




7.819 


3.544 


2.454 


1.762 


1.076 


0.591 


0.340 


0.148 


0.111 


0.083 


0.062 


0.053 


/r, % 


23.9 


18.5 


18.3 


20.0 


16.2 


17.5 


19.3 


19.9 


20.4 


21.2 


22.3 


23.4 




5.776 


3.243 


2.516 


1.883 


1.184 


0.846 


0.480 


0.223 


0.166 


0.100 


0.085 


0.081 




8.201 


4.221 


3.151 


2.320 


1.426 


1.023 


0.587 


0.272 


0.203 


0.130 


0.111 


0.105 


/s, % 


29.6 


23.2 


20.1 


18.8 


17.0 


17.3 


18.1 


17.9 


18.0 


23.3 


23.0 


22.5 




Figure 3: The cross sections of Coulomb deexcitation cr^° multiplied by Ecm for {fip)n + H 
collisions calculated in the quantum-mechanical close-coupling method in comparison with the 
results of Bracci and Fiorentini [11] (B&F). 



influence of the open channels with An > 1 which increases for the higher n states and can't be 
taken into account in the two-state approaches [11-13]. 

The distribution over the final states n' is completely different from the semiclassical results 
[11] as illustrated in Fig. 4 and Table 3. The present calculations predict that An = 1 transitions 
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dominate in agreement with the semiclassical, AAA and CTMC results however the transitions 
with An > 1 are strongly enhanced as compared to the results [11-13]. To determine the fraction 
of the transitions with An > 1 in the total CD cross section, the values 



fn 



<^n ~ '^n,n-l 



cr, 



CD 



100%. 



are also presented. 

The calculated CD cross sections an,n-i{E), a^^{E) and fn{E) for {fip)n + H colhsions with 
n = 3 -h 8 and energies from 0.01 up to 20 eV are presented in Table 3. The observable variations 
of the CD cross sections with E can be related to the opening of the additional channels with 
An > 1. The violation of the known in literature [11-13] n-dependence of the CD cross sections is 
illustrated in Figs. 2 and 3 and in more detail in Table 3. For n = 4 6 we observe the essential 
variations in the CD cross sections in contrast to the two-state calculations [11-13]. In Fig. 3 
this effect can be clearer because of the different energy dependence of the CD cross sections for 
various n. In our opinion the traditional n dependence {cr^^ ~ n'' with 7 > 2) followed from 
the semiclassical picture can be reached at the relative energies comparable or much more than 



3 



0.001 




T 

0° 0.014— 



0.01 



0.001 



Figure 4: The Coulomb deexcitation cross sections an,n' in a.u. for {fip)n+H collisions a.t E = 1 eV 
calculated in the quantum-mechanical close-coupling method in comparison with results of Bracci 
and Fiorentini [11]. 



The transitions with An = 1 are most likely but the An > 1 transitions make up a substantial 
fraction of the total CD cross section (16% - 37%) for n > 4 in contrast to the two-state approaches 
[11-13]. The energy dependence of this fraction is completely different for various n. In particular, 
this behaviour leads to the unexpected dependences of the cross sections on the principle quantum 
number n and coUisional energy -Ecm- 

It is possible that the problem of the high energy fraction of the kinetic energy distribution 
in the muonic hydrogen is related to the unproper behaviour of the CD cross sections taken in 
the extended standard cascade model (ESCM) [17] for n = 3 -j- 7 and based on the semiclas- 
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sical model [11]. The dependences of the CD cross sections on n and E used in ESCM are in 
disagreement with the results of our more elaborate close-coupling consideration. 




E ,eV 

cm 

Figure 6: The cross sections of Coulomb deexcitation cr„,„_i in a.u. for {^d)n + H collisions 
calculated in the close-coupling approach (in two-level approximation) (CC) in comparison with 
results of Ponomarev and Solovyov [13] (P&S). 

In Fig. 6 we compare the present cross sections of CD an,n~i{E) for + H collisions with 

the results of AAA [13] for n = 3,5. As in case of {^p)n + H [11] the transition 3 2 is strongly 
suppressed (almost two order) in comparison with the close-coupling calculations. At the same 
time for the transition 5 — > 4 the results of AAA are only suppressed in four times as compared 
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with our results. The reason of this discrepancy is not clear at present. 



4 Conclusion 

The unified treatment of the elastic scattering, Stark transitions and Coulomb deexcitation is 
presented in ab initio quantum-mechanical approach. The main features of the Coulomb deex- 
citation in the collision of the exotic hydrogen atom in excited states with hydrogen atom have 
been investigated in detail using the close-coupling method. The new results for the n and E 
dependences of the CD cross sections are obtained for n = 3 -i- 8 and relative energies up to 20 
eV relevant to the kinetics. The calculated cross sections do not agree with the long-standing 
traditional beliefs about a^^{E) ~ rC /E dependence of CD cross sections which is based on 
semiclassical approximation. The new important results are also obtained for the fraction of the 
transitions with An > 1. It is shown that the contribution of these transitions is more essential 
than it is assumed earlier and reaches up to ~37% for n = 6. We suppose that our results pro- 
vide a reliable theoretical input for the further kinetics calculations. More detailed investigations 
concerning the partial transitions and other exotic particles will be discussed in future publications. 
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